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ABSTRACT 

We present a critical assessment of commonly used pre-main-sequence isochrones by 
comparing their predictions to a set of well-calibrated colour-magnitude diagrams 
of the Pleiades in the wavelength range 0.4 to 2.5 /am. Our analysis shows that for 
temperatures less than 4000 K the models systematically overestimate the flux by a 
factor two at 0.5 /im, though this decreases with wavelength, becoming negligible at 
2.2 /mi. In optical colours this will result in the ages for stars younger than lOMyr 
being underestimated by factors between two and three. 

We show that using observations of standard stars to transform the data into 
a standard system can introduce significant errors in the positioning of pre-main- 
sequences in colour-magnitude diagrams. Therefore we have compared the models to 
the data in the natural photometric system in which the observations were taken. Thus 
we have constructed and tested a model of the system responses for the Wide-Field 
Camera on the Isaac Newton Telescope. 

As a benchmark test for the development of pre-main-sequence models we provide 
both our system responses and the Pleiades sequence. 

Key words: stars: evolution - stars: formation - stars: pre-main-sequence - stars: 
fundamental parameters - techniques: photometric - open clusters and associations: 
general - Hertzsprung-Russell and colour-magnitude diagrams 



1 INTRODUCTION 

Our understanding of time-dependent physical processes in 
young pre-main-sequence (pre-MS) and main-sequence (MS) 
populations is fundamentally limited by poorly constrained 
timescales, the derivation of which requires the determina- 
tion of ages for pre-MS stellar populations. These ages are 
generally derived via the comparison of photometric obser- 
vations of pre-MS stars in young open clusters with theoreti- 
cal models (or isochrones) using colour-magnitude diagrams 
(CMDs). 

Many previous studies deriving stellar parameters from 
observations of young clusters in this way have shown 
large discrepancies between the photometric data and the 
model predictions. It has been well documented that pre- 
MS isochrones generally do not fit the entire dataset in 
CMDs; e.g. whilst tracing the sequence defined by the higher 
mass members they will not simultaneously follow the se- 
quence in the lower mass regime (see for instance Hartmann 
2003; Stauffer et al. 2007). In addition, pre-MS isochrones 
from different evolutionary models give ages that differ by 
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a factor two for the same cluster (e.g. Dahm 2005; Mayne 
et al. 2007). Furthermore, the same dataset compared to 
isochrones in different colours can lead to different age de- 
terminations (e.g. Naylor et al. 2002). 

The situation has become even more uncertain with the 
realisation that ages derived from the MS are systematically 
older by a factor 1.5 to 2 than ages derived from the pre-MS 
(e.g. Naylor 2009). To make further progress we require a 
precise test of the pre-MS models to determine whether the 
ages could be in error by such a large factor, and if any parts 
of the pre-MS or any colour combinations are reliable age 
indicators. Our longer term aim is that if we have a clear 
understanding of the limitations of pre-MS ages we may be 
able to create a reliable pre-MS age scale. 

In this work, therefore, we set up a benchmark test for 
pre-MS isochrones. Most obviously that benchmark must 
be a cluster that contains a significant number of pre-MS 
objects, as well as a populated MS to show that the models 
at least fit these more evolved stars. In addition it should 
have distance and age determinations which are independent 
of the CMD, and preferably a low extinction (which might 
otherwise complicate transformations) . Given that most pre- 
MS isochrones are only available for solar metallicity, the 
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cluster must be at least close to solar composition. Finally, 
it should be nearby so that the lowest possible masses can be 
accessed. The Pleiades is the only cluster which fits all these 
criteria, and Stauffer et al. (2007) have already undertaken 
a comparison in V and I c , which shows that at least the 
V-band luminosity is not matched by two models available 
at that time. 

In this paper we return to the Pleiades with the aim of 
assessing which of the most recent models best match the 
data. This then allows us to quantify the remaining mis- 
match so that we can make a realistic assessment of its im- 
pact on pre-MS ages. We trace the mismatch as a function 
of wavelength using observations in a set of filters which are 
contiguous in wavelength and cover 0.4 to 2.5 /xm. The data 
reach cooler temperatures than the Stauffer et al. (2007) 
study, and are free of the uncertainties introduced by trans- 
forming heterogeneous datasets into a standard photometric 
system. 

In Section 2 we describe the set of interior and atmo- 
spheric theoretical models we have studied and examine sys- 
tematic differences between the different formalisms. In Sec- 
tion 3 we describe a set of observations, of both standard 
stars and the Pleiades, which we use to test the models. 
In Section 4 we compare the models to the Pleiades dataset 
using CMDs. Section 5 describes the comparison of the mod- 
els to a sample of MS binaries with dynamically constrained 
masses. In Section 6 we quantify the mismatch between the 
models and the data as a function of wavelength. Our con- 
clusions are detailed in Section 7. 



2 THE MODELS 

To derive fundamental parameters, such as age and distance, 
from fitting CMDs we require model isochrones, which must 
also be transformed into the requisite photometric system. 
Isochrones are constructed from stellar interior models that 
predict Lboi, ^eff and surface gravity (\ogg). Colours are 
then calculated through a colour- T e ff relation and magni- 
tudes via bolometric corrections to Lboi- Both relations can 
be derived by folding atmospheric model flux distributions 
through the appropriate photometric filter responses. The 
colours and magnitudes must then be calibrated to a stan- 
dard scale. Therefore, for this study we have adopted a set of 
pre-MS interior models and also a collection of atmospheric 
models with which to derive the required relations. 

2.1 Pre- main-sequence interior models 

We tested those pre-MS interior models that were publicly 
available, covered a wide range in stellar mass and had been 
designed specifically for the investigation of pre-MS evo- 
lution. The sets of interior models studied here are those 
of Baraffe et al. (1998), Siess, Dufour, & Forestini (2000), 
D'Antona & Mazzitelli (1997) and Dotter et al. (2008) (here- 
after BCAH98, SDF00, DAM97 and DCJ08 respectively 1 ). 
Note that the recent Pisa models (Tognelli, Prada Moroni, & 



1 Whilst the DCJ08 models are not specifically designed for 
studying pre-MS evolution, these mass tracks cover the entire 
pre-MS phase over a significant mass range. 
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Figure 1. Variation between the pre-MS evolutionary tracks for 
masses of 0.1,0.5, and 1.0 Mq (top panel) and 4.6 Gyr isochrones 
(bottom panel) for the following models: BCAH98 a = 1.9 
(red, continuous), BCAH98 a = 1.0 (red, dashed), SDF00 
(blue), DAM97 (green) and DCJ08 (cyan). Note that for masses 
< 0.6 Mq both computations of the BCAH98 models are iden- 
tical, this is due to the insensitivity to the value of the mixing 
length parameter in this mass regime. 

Degl'Innocenti 2011) do not extend to the age of the Pleiades 
(see Section 4.1) and therefore are not used in this study. 

Fig. 1 shows the variation between the different pre-MS 
evolutionary tracks and the 4.6 Gyr isochrone for a range of 
masses. There are systematic differences in both the mass 
tracks, especially at young ages, and location of the MS, 
particularly for low-mass stars. These differences stem from 
variations in the treatment of various physical processes in 
addition to the values of adopted parameters. Most notable 
are the treatment of convection, the opacity sources and the 
treatment of the stellar interior/atmosphere boundary (see 
Hillenbrand & White 2004). 

The BCAH98 models have been computed for two dif- 
ferent mixing length parameters a — 1.0 (general mixing 
length parameter) and a = 1.9 (solar calibrated value), 
where for stellar masses below 0.6 M© these models are iden- 
tical as the dependence on the mixing length parameter is 
negligible in this regime (Baraffe et al. 2002). Further note 
that the publicly available SDF00 models have a slightly 
modified solar abundance and mixing length parameter from 
the solar calibrated model (cf. Z — 0.02 and a = 1.5 with 
Z = 0.0249 and a = 1.605 in Siess et al. 2000). 



2.2 Atmospheric models 

We required atmospheric model flux distributions to calcu- 
late a theoretical colour- T e ff relation and bolometric cor- 
rections. We produced a spectral library consisting of the 
PHOENIX BT-Settl model atmospheres for 400 < T e ff < 
7800 K and the ATLAS9 models with newly updated opac- 
ity distribution functions (ODFnew) for 8000 < T e ff < 
50 000 K. Despite the differences in microphysics between 
the two sets of atmospheric models, we found that at the 
transitional T e ff = 8000 K derived colours from both sets of 
models agree to within 0.02 mag in all colours. 
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2.2.1 The ATL AS9 /ODFnew models 

The ATLAS 9 /ODFnew atmospheric models are those of 
Castelli & Kurucz (2004) 2 . These models are based on the 
assumption of steady-state plane-parallel layers under the 
influence of local thermodynamic equilibrium and computed 
using the solar abundances of Grevesse & Sauval (1998). 
Line blanketing effects are computed statistically via the 
use of the opacity distribution functions which average the 
contribution from the various atomic and molecular species 
as described in Kurucz (1979). A pure mixing length theory 
(Bohm-Vitense 1958) is used in the computation of these 
models and we adopted those termed "no-overshoot" on the 
basis that these better recreate the observed spectra of stars 
with T e ff greater than the Sun (Castelli, Gratton, & Kurucz 
1997). The mixing length parameter is set to a — 1.25 with 
a microturbulent velocity £ = 2kms -1 . 

2.2.2 The BT-Settl models 

The BT-Settl atmospheric models are those described in 
Allard, Homeier, & Freytag (201 1) 3 . These models have 
been computed using an updated version of the phoenix 
atmospheric code (cf. Allard et al. 2001) to include the up- 
dated Barber et al. (2006) BT2 H 2 line list, the revised 
solar abundances of Asplund et al. (2009) and a sophisti- 
cated cloud model that accounts for the settling of dust 
grains (Allard et al. 2003). The effects of turbulent mix- 
ing in the atmosphere are calculated by interpolating from 
2-D and 3-D radiation hydrodynamic models (Ludwig, Al- 
lard, & Hauschildt 2006), with the treatment of dust for 
cooler models (T e ff < 2600 K) based on the dust formation 
models of Freytag et al. (2010). The effects of opacity are 
sampled directly using a library of over 700 million lines 
including atomic and molecular species along the spectrum. 
The BT-Settl models are computed under the assumption of 
plane-parallel radiative transfer, where convection is treated 
using the standard mixing length theory. The mixing length 
parameter is defined to be a — 2.0 with a microturbulent 
velocity of £ = 2kms _1 . 

3 THE DATA 
3.1 Observations 

Our observations are a subset of a large survey of pre-MS 
clusters (the remainder of which will be published in an up- 
coming paper; Bell et al. in prep.) and were obtained using 
the 2.5-m Isaac Newton Telescope (INT) on La Palma. The 
survey was split over two runs, using the same instrumenta- 
tion and filter set on both occasions, namely the four EEV 
2k x 4k CCD Wide-Field Camera (WFC) with a 34' x 34' 
fleld-of-view and the (UgriZ) WFC filter set. A combination 
of long and short exposures were used to ensure that any 
bright stars saturated in the longer exposures could be mea- 
sured in the shorter exposures. Whilst the (gri) WFC filters 
are based on the Sloan Digital Sky Survey (SDSS) design 
there are no similar u and z filters for the WFC and so 

2 ftp : //ftp . stsci . edu/cdbs/grid/ck04models 

3 http: //www. phoenix. ens-lyon.fr/Grids/BT-Settl/ 
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Figure 2. Two square degree three-colour mosaic image of the 
Pleiades from the ESO Digital Sky Survey with the five INT-WFC 
fields-of-view overlaid. Note that the image is slightly displaced 
from the central coordinates of a = 03 h 47 m 24 s , 5 = +24° 07' 00" 
to illustrate the region chosen to identify low-mass members. 

the Royal Greenwich Observatory (RGO) U WFC and Z WFC 
filters were used. The first set of observations were taken 
during the 21 st - 26 th Oct. 2007 with the second run during 
the 13 th - 19 th Sept. 2008. 

Five overlapping fields within the Pleiades open cluster 
were observed over the course of the two runs (see Fig. 2 
and Table 1). In addition, four fields in Stripe 82 of the 
SDSS were used as standard fields and routinely observed 
throughout both runs. 

3.2 Data reduction 

Each image was debiased and flat fielded using a median 
stacked bias and flat field frames. Both the i WFC - and Z WFC - 
band images were def ringed using library fringe frames. The 
data were analysed using the cluster package as described 
in Naylor et al. (2002), Burningham et al. (2003) and Jef- 
fries et al. (2004). The i WFC -band images for a particular 
field were combined (after calculating spatial transforma- 
tions between separate frames) with the resulting deep im- 
age used for object identification and detection. An object 
list was created and optimal photometry (Naylor 1998) was 
then performed providing a flux for each object at a given 
position. A spatially varying profile correction was calcu- 
lated from the bright unsaturated stars in each frame so 
that the photometry could later be calibrated using stan- 
dard star observations. Although the i WFC -band images were 
combined, optimal photometry was carried out on each indi- 
vidual frame, so as to ensure that the good signal-to-noise in 
one frame was not outweighed by a poorer signal-to-noise in 
another frame upon combining the images. The individual 
profile corrected measurements of each frame were then com- 
bined by weighting each measurement in accordance with its 
signal-to-noise ratio and accounting for differences in air- 
mass between separate frames. An additional statistical un- 
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Table 1. The central coordinates for each neld-of-view and ex- 
posure times in the INT-WFC bandpasses for the observations of 
the Pleiades. 



Field 


RA Dec. 


Filter 


Exposure time (s) 




(J2000.0) 




x 1 unless stated 


Pleiades 


03 h 46 m 19.2 s 


TJ 

WFC 


1, 10 


Field A 


+23° 38' 02.4" 


^WFC 


1, 10, 100, 500(x 2) 






r WFC 


1, 7, 50, 250 






Vfc 


1, 10, 100(x 2) 








1, 7, 50, 250 


Pleiades 


03 h 45 m 00.0 s 


TJ 

WFC 


1, 10 


Field B 


+23° 44' 45.6" 


S'WFC 


1, 10, 100, 500(x 2) 






r WFC 


1, 7, 50, 250 






Vfc 


1, 10, 100(x 2) 








1, 7, 50, 250 


Pleiades 


03 h 43 m 44.4 s 


TJ 

WFC 


1, 10 


Field C 


+23° 58' 08.4" 




1, 10, 100, 500(x 2) 






r WFC 


1, 7, 50, 250 






Vfc 


1, 10, 100(x 2) 






\fc 


1, 7, 50, 250 


Pleiades 


03 h 42 m 36.0 s 


^WFC 


1, 10 


Field D 


+23° 58' 08.4" 


S'WFC 


1, 10, 100, 500(x 2) 






r WFC 


1, 7, 50, 250 






Vfc 


1, 10, 100(x 2) 






\fc 


1, 7, 50, 250 


Pleiades 


03 h 43 m 44.4 s 


^WFC 


1, 10 


Field E 


+24° 24' 28.8" 


S'WFC 


1, 10, 100, 500(x 2) 






r WFC 


1, 7, 50, 250 






Vfc 


1, 10, 100(x 2) 






^WFC 


1, 7, 50, 250 



cess ensures greater consistency between fields in the final 
catalogue. The resulting zero-point shift is an indicator of 
the internal consistency of the photometry, as well as the 
accuracy with which the profile corrections was performed, 
and suggests an accuracy of c± 1 per cent. Our full Pleiades 
photometric catalogue is given in Table 2. 

3.3 Zero-point stability 

Standard field observations were taken on five separate 
evenings. On only one of the five nights were conditions 
photometric throughout, the other four nights were af- 
fected either at the beginning or the end by the pres- 
ence of cloud. We therefore calculated a zero point for 
each g WFC -band observation by comparing the instrumen- 
tal magnitudes to those in the SDSS catalogue, and cor- 
recting for extinction using a mean g WFC -band extinction 
coefficient (k Q = 0.19) for La Palma. We could then 
identify when the zero-point became unstable and define 
half nights of bona fide photometric conditions, which were 
then used for the photometric calibration procedure (the 
maximum deviation over the course of a single half night 
was 0.05 mag). Due to a paucity of blue stars in our stan- 
dards catalogue, we opted to use the colour range which 
was most populated (0.6 < (g — i) WFC < 1.0) and define 
our zero-point at a median colour of (g — i) WFC =0.8. This 
colour range cut in (g — i) WFC defines the necessary cuts in 
the other INT-WFC colours (1.05 < (U - g) WFC < 1.96, 
0.14 < (r - i) WFC < 0.31 and 0.02 <(i- Z) WFC < 0.18 with 
median values of (U — g) WFC = 1.505, (r — i) WFC = 0.225 
and (i — Z) WFC = 0.1). The seeing generally varied between 
1 and 2 arcsec in the temporal regions defined above. 



certainty was added to each measurement at this stage to 
ensure that the distribution of chi-squared (% 2 ) resulting 
from combining the measurements versus signal-to-noise was 
independent of signal-to-noise and of order unity (Naylor 
et al. 2002). The additional statistical uncertainty adopted 
here ranged from 0.01 — 0.03 mag and reflects uncertainties 
in the profile correction. Stars with a reduced x 2 > 10 are 
generally nagged as variable (see Burningham et al. 2003). 
For some fields-of-view, we have combined photometric ob- 
servations at two different epochs and due to pre-MS vari- 
ability over such timescales, there are a number of sources 
that have thus been nagged. For sources which displayed 
night- to- night variations, we have replaced the photometry 
with colours and magnitudes from a single night. The Two- 
Micron All-Sky Survey (2MASS, Cutri et al. 2003) was used 
to provide an astrometric solution for each frame. The rms 
residual of the 6-coefficient fit were approximately 0.2 arcsec. 

To create a final optical catalogue of the Pleiades based 
on the five overlapping fields, it was necessary to combine the 
individual fields-of-view. For this process, we used the nor- 
malisation procedure as described in Jeffries et al. (2004). A 
running catalogue was created that contained all the derived 
magnitudes and colours for each object identified in all fields. 
We then calculated the mean magnitude and colour differ- 
ences of stars in common between two overlapping fields, 
and adjusted the zero-points for each field to minimise these 
differences (in effect allowing for small variations in the zero- 
point between fields observed on different nights) . This pro- 



3.4 Transforming the data into the standard 
system 

Traditional photometric calibration transforms instrumen- 
tal magnitudes into magnitudes on a standard system. Typ- 
ically colours and magnitudes are transformed using linear 
functions of the form 

#stand =1pg(g- i)inst ~ k g X + Z g , (1) 
(g - ^)stand = 1pgi(g ~ «)inst ~ kgiX + Z g i, (2) 

where x is the airmass, k the extinction coefficients, z the 
zero-points and ip the colour terms (k, z and ip are deter- 
mined from observations of standard stars). Fig. 3 shows 
the g-band magnitude and g — i colour differences between 
the standard SDSS and the natural INT-WFC photometric 
systems (see Appendices A and B for details). It is immedi- 
ately obvious that there is no linear colour-dependent trans- 
formation that can be applied to INT-WFC observations of 
MS stars to calibrate them into the standard SDSS photo- 
metric system. Even if non-linear MS transformations were 
created, it is clear that they could not be used to trans- 
form observations of red ((g — i)sDSS > 1.8) pre-MS stars 
into the standard system. This is due to differences in the 
spectra between MS and pre-MS stars of the same colour. 
These differences can result in both g-band magnitude and 
g — i colour errors of order 0.1 — 0.15 mag at a colour of 
(g — i) sr> ss — culminating in stars that ultimately occupy 
the wrong position in CMD space. Even larger errors can be 
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Table 2. A sample of the full Pleiades photometric catalogue with colours and magnitudes in the natural INT-WFC photometric system. 
Due to space restrictions, we only show the g WFC and (g — i) WFC colours and magnitudes as a representation of its content. The full 
photometric catalogue (available online) also includes photometry in the (U — g) WFC , (r-i) WFC , and (i — Z) WFC colours and the U WFC , 
r WFC , i WFC , and Z WFC magnitudes. Columns list unique identifiers for each star in the catalogue: field and CCD number (integer and 
decimal), ID, RA and Dec. (J2000.0), CCD pixel coordinates of the star, and for each of g WFC and (g — i) WFC there is a magnitude, 
an uncertainty on the magnitude and a flag (OO represents a "clean detection" ; see Burningham et al. 2003 for a full description of the 
flags.) 



Field 


ID 


RA (J2000.0) 


Dec. (J2000.0) 


X 


y 


^WFC 


(Ja 

g WFC 


Flag 


(9 ^)wfc 




Flag 


21.04 


8 


03 46 34.136 


+23 37 26.25 


1204.056 


1439.980 


8.236 


0.010 


OO 


-0.412 


0.014 


OL 


22.02 


16 


03 43 41.521 


+23 38 56.30 


998.324 


2035.800 


8.523 


0.010 


OS 


-2.653 


0.014 


ss 




(g-i)s 



Figure 3. Calculated transformations between the SDSS and 
INT-WFC photometric systems. The transformations have been 
calculated for reddened (nominal E(B — V) = 1 mag; red) and 
unreddened (blue) ZAMS (crosses), lOMyr (asterisks) and 1 Myr 
pre-MS stars (open circles). Upper panel: Ag versus (g — i) SJDSS - 
Lower panel: A(g — i) versus (g — i) sr>ss - No single (even non- 
linear) transformation will calibrate observations of both red ({g— 



> 1.8) MS and pre-MS stars into a standard system, and 



if either are reddened the situation becomes worse (c 
at a colour (g — i) SDSS — 3). 
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Figure 4. Magnitude difference in the t/ WFC -band between the 
WFC ca i c and WFC b s photometric catalogues as a function of 
(U — fiOwFc-caic- Asterisks represent the median value of all points 
in a given bin of size 0.25 mag in colour. The error bars on these 
symbols are the standard error about the median value. Open 
circles represent individual data points and are plotted where the 
number of points is insufficient to calculate the median value in a 
given bin (defined as five). The error bars of these symbols are the 
individual uncertainties associated with that point. The dashed 
lines represents the ±0.02 mag level with respect to zero. 



caused by the difference between the spectrum of a reddened 
and an unreddened star of the same apparent colour (see 
Fig. 3). From the calculated SDF00 zero-age main-sequence 
(ZAMS; see Appendix B) the differential behaviour between 
MS and pre-MS transformations occurs at an approximate 
spectral type of between K6 and M0. Mayne et al. (2012) 
document a similar, though more extreme situation which 
can result from using MS transformations for brown-dwarfs. 

A deviation of up to 0.15 mag may help explain, to some 
extent, why pre-MS isochrones do not simultaneously fit 
both the higher and lower mass members within a given 
cluster (Hartmann 2003; Stauffer et al. 2007; Jeffries et al. 
2009), although other factors such as photometric variabil- 
ity, extinction, unresolved binaries, accretion luminosity and 
poorly constrained opacity line lists in the atmospheric mod- 
els may also contribute to this (Hartmann 2001; Burning- 
ham et al. 2005). Furthermore, this could also have severe 
implications for cluster age and mass function determina- 
tions. Allowing the pre-MS stars in a CMD to "float" by 
0.15 mag in both colour and magnitude can result in an 



age difference of a factor of two through isochrone fitting. 
The conversion from observable magnitudes and colours into 
mass estimates are highly age and model dependent. Hence 
an incorrrect age would result in an erroneous mass func- 
tion estimate thus making meaningful comparison between 
cluster mass functions difficult at best. 



3.5 Calibrating the natural system 

If the observations cannot be transformed into a standard 
photometric system, then the photometry must remain in 
the natural photometric system and the theoretical mod- 
els transformed into this system. To achieve this we need 
a reliable model for the throughput for each bandpass as a 
function of wavelength. We emphasise that this is not sim- 
ply a matter of using the filter responses, or even multiply- 
ing the CCD quantum efficiency by the filter response, but 
involves allowing for every element from the Earth's atmo- 
sphere through to the CCD to create a system response. 
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Figure 5. Same as Fig. 4 but for the difference in the g WFC -band 
as a function of (g - *) WFC;calc . 
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Figure 8. Same as Fig. 4 but for the difference in the Z WFC -band 
as a function of (i - ^) WFC;calc - 
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Figure 6. Same as Fig. 4 but for the difference in the r WFC -band 
as a function of (r - ^) WFC;calc - 




1 2 

calc 

Figure 7. Same as Fig. 4 but for the difference in the i WFC -band 
as a function of (g - *) WFC;calc . 



We present our best estimates for the system responses in 
Appendix A. 

Having obtained the system responses it is crucial that 
we obtain an estimate of how precise these are. In prin- 
ciple we could achieve this by folding the observed fluxes 
of spectrophotometric standards through our estimated sys- 
tem responses and seeing how closely the observed and pre- 
dicted magnitudes match. In practice such standards have 
not been measured in sufficient numbers with the range of 
colours and magnitudes to prove useful for our purposes. 
Instead we took a sample of MS stars, used their observed 
SDSS colours to establish their effective temperatures, and 
then folded model atmospheres of the appropriate temper- 
ature through our estimated system responses to find their 
magnitudes in the INT-WFC system. 

We defined our sample of MS stars using a series of 
colour-colour diagrams to isolate the observed MS in the 
Stripe 82 fields and trim any obvious outliers. To obtain 
colours and magnitudes in the INT-WFC system we then 
transformed the standard SDSS colours and magnitudes 
into the INT-WFC photometric system (hereafter WFC ca ic). 
These transforms are calculated by folding model atmo- 
spheres through both the SDSS and our estimated INT- 
WFC responses and are given in Table 3. The details of these 
calculations are given in Appendix B, and the transformed 
standards in the INT-WFC colours are given in Table 4. 

With the INT-WFC standards established, we can eval- 
uate how well constrained the calculated system responses 
of the INT-WFC are. Calculating zero-points in the colour 
ranges specified in Section 3.3 and applying these zero-points 
and extinction coefficients only (no colour terms) to all stan- 
dard field observations we created a set of observed stan- 
dard star colours and magnitudes in the INT-WFC system. 
To create a single catalogue for each of the four standard 
fields we performed the normalisation process as described 
in Section 3.2 and merged all four fields to create a catalogue 
of observed colours and magnitudes of Stripe 82 standards 
in the INT-WFC photometric system (hereafter WFC G bs)- 
After the normalisation procedure, we examined, for each 
star, the magnitude shift between each observation and the 
mean observed magnitude. This is an indicator of the inter- 
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Table 3. The calculated transformations between the SDSS and INT-WFC photometric systems for unreddened MS stars where A 
represents (INT-WFC — SDSS). The full table is available online, a sample is shown here as a representation of its content. 



T e fi 




Ag 


A(g - i) 


(9 OsDSS 


A(u-g) 


( U SOsDSS 


A(r - i) 


( r ^SDSS 


A(i-z) 


(* Z )sDSS 


2775 


5.25 


-0.29300 


-0.04765 


3.83961 


0.01264 


4.12961 


0.24688 


2.20001 


-0.17476 


1.22123 


3290 


5.11 


-0.23094 


-0.10391 


2.68386 


-0.02563 


3.06419 


0.12898 


1.20809 


-0.08704 


0.66867 


3528 


4.98 


-0.21237 


-0.11299 


2.35355 


-0.03358 


2.87477 


0.09764 


0.94887 


-0.06707 


0.53108 


3683 


4.92 


-0.20428 


-0.12031 


2.19344 


-0.03595 


2.83717 


0.07940 


0.81597 


-0.05550 


0.45615 


3829 


4.87 


-0.19761 


-0.12806 


2.05653 


-0.03751 


2.81937 


0.06188 


0.70141 


-0.04429 


0.38859 



Table 4. A sample of the transformed Stripe 82 standard star catalogue with colours and magnitudes in the natural INT-WFC photo- 
metric system. The full photometric (available online) also includes photometry in the (U — g) WFC , (r — i) WFC , and (i — Z) WFC colours. 
Columns list unique identifiers for each star in the catalogue (see Table 2) for details. 



Field 


ID 


RA (J2000.0) 


Dec. (J2000.0) 


X 


y 


P-WFC 


^WFC 


Flag 


(9 ^)wfc 


<j( a -j\ Flag 

\9 ^WFC 


0.00 
0.00 


67737 
67740 


20 47 07.346 
20 47 09.412 


-01 05 52.53 
-01 04 17.09 


0.000 
0.000 


0.000 
0.000 


17.158 
14.721 


0.006 
0.006 


OO 
OO 


0.585 
0.482 


0.008 OO 
0.008 OO 



nal precision of our photometry and we found that g rms = 
0.011 mag, (g - i) rms = 0.014 mag, (U - #) rm s = 0.020 mag, 
(r — i) rms = 0.011 mag and (i — Z) rms — 0.010 mag. We at- 
tribute the slightly poorer precision in U WFC to a lack of 
colour dependent atmospheric extinction terms. 

The differences WFC ca ic — WFCobs are shown in 
Figs. 4-8. Using bin sizes of 0.25 mag in colour, and only 
data with uncertainties of less than 0.03 mag we calculated 
the median value in each colour bin as an indicator of the 
level of agreement between the WFC ca ic and WFCobs cata- 
logues in a given colour range. A paucity of data at extreme 
colours forced us to relax our restrictions on the uncertainty 
in the data points to less than 0.1 mag and where the num- 
ber of stars was insufficient to calculate the median value in 
a given colour bin (taken as five) , the individual data points 
were plotted. These plots show that we reproduce the zero 
point to within 0.01 mag at the colours which we denned 
to calculate the zero-points in Section 3.3. The remaining 
residuals lie within the ±0.02 mag level across the entire 
colour range except for a small region in the £/ WFC -band 
and colours redder than (g — i) WFC — 2 in the i WFC -band. 

The most obvious causes of any disagreement between 
WFCcaic and WFCobs are that our estimated system re- 
sponses are incorrect, and that the atmospheric models 
are incorrect. In an attempt to obtain some idea of the 
likely magnitude of the latter effect we tried re-calculating 
the transforms using different model atmospheres (the 
PHOENix/GAIA models of Brott & Hauschildt 2005) and 
the observed spectral library of Pickles (1998). We found 
these changed the transformations by up to 0.02 mag, nicely 
explaining the majority of the residuals. It remains unclear 
whether the remaining residuals are due to errors in the 
atmospheric fluxes, or poorly modelled U WFC - and i WFC - 
band system responses. As a result, we retain the system 
responses in their current calculated form for the remainder 
of this study, and note the fact there are uncertainties at the 
levels shown in Figs. 4-8. 



4 COMPARING THE MODELS AND THE 
DATA - THE PLEIADES 

Having characterised the natural INT-WFC photometric 
system and calculated the corresponding system responses, 
we can now test the pre-MS interior and stellar atmo- 
spheric models by comparing them to our Pleiades data. 
We obtained a list of approximately 180 members, by cross- 
correlating our full catalogue with the combined member- 
ship catalogues of Stauffer et al. (2007) and Lodieu, Deacon, 
& Hambly (2012). This membership list is given in Table 5. 

4.1 Model parameters 

We used a distance modulus dm = 5. 63 ±0.02 which is based 
on the trigonometric parallax of Soderblom et al. (2005). We 
chose to rely on the lithium depletion age, since although it 
is based on the predictions of stellar interior models, Jeffries 
& Oliveira (2005) show there is a very high level of agree- 
ment between the various models. Furthermore we used the 
age based on the boundary in K s (130 Myr) of Barrado y 
Navascues, Stauffer, & Jayawardhana (2004), since as we 
shall show later, K s is the best predictor of luminosity avail- 
able to us. We used a reddening of E(B — V) = 0.04 based 
on the mean extinction Ay — 0.12 (see Stauffer et al. 1998) 
and Ry — 3.2. Hence, to transform the theoretical models 
into the observable plane, we first reddened the atmospheric 
models by a nominal E(B — V) = 0.04 using the reddening 
law of Cardelli, Clayton, & Mathis (1989) and then calcu- 
lated bolometric corrections using Eqn. B2. 

4.2 Discussion 

Fig. 9 shows four CMDs of the Pleiades with the pre-MS 
isochrones of BCAH98, SDF00, DAM97 and DCJ08 over- 
laid. The break between MS and pre-MS stars in the Pleiades 
occurs at (g — i) WFC — 2, and is marked by a paucity of 
stars in all CMDs. It is clear that the isochrones all fol- 
low the Pleiades MS reasonably well, though as we shall see 
in Section 6 this is in part is due to compensating errors 
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Table 5. A sample of the catalogue of colours and magnitudes in the natural INT-WFC photometric system for Pleiades members. The 
columns and content are in the same format as Table 2. The full photometric catalogue is available online. 



Field 


ID 


RA (J2000.0) 


Dec. (J2000.0) 


X 


y 


S\VFC 


%FC 


Flag 


(9 ^)wfc 


a (g-i) WFC 


Flag 


24.02 


45 


03 41 03.006 


+23 43 21.42 


426.717 


3696.221 


17.653 


0.006 


oo 


2.553 


0.006 


OO 


24.02 


118 


03 41 05.230 


+23 50 14.82 


526.182 


2444.212 


20.400 


0.015 


oo 


3.228 


0.017 


00 




1 2 3 4 0.2 0.4 0.6 0.8 1.0 1.2 

Figure 9. Optical/near-infrared CMDs of Pleiades members. The 130 Myr pre-MS isochrones of BCAH98 a = 1.9 (red, continuous), 
BCAH98 a = 1.0 (red, dashed), SDF00 (blue), DAM97 (green) and DCJ08 (cyan) are overlaid, adopting a distance modulus dm = 5.63 
and a reddening E(B — V) = 0.04. The isochrones have been transformed into the observable plane using bolometric corrections derived 
by folding the atmospheric model flux distributions (see Section 2.2) through the calculated INT-WFC system responses. The black line 
in each panel represents the spline fit (by eye) to the Pleiades single-star sequence. The open circles represent the photometric data, with 
the associated uncertainties in colour and magnitude shown as the bars. Top left: K s ,g WFC — K s CMD. Top right: g WFC , (g — i) WFC 
CMD. Bottom left: r WFC , (i - Z) WFC CMD. Bottom right: K s , J — K s CMD. 



in each individual band. None of the models, however, trace 
the observed pre-MS locus, with each predicting colours that 
are too blue for a given stellar mass. This issue has already 
been highlighted by Stauffer et al. (2007) (see their Fig. 13), 
though from their data one could only be certain there was 
a problem in the V-band, whereas our data show that all 
optical bands are affected. 

This is part of a broader picture where it is known 
that MS stars are similarly affected, at least in the l/-band 
(see Section 5). It is well known that there are sources 
of opacity missing from the atmospheric models for cool 
stars (T e ff < 3700 K) in the optical regime of the spectrum 



(Leggett et al. 1996; Baraffe et al. 1998; Alvarez & Plez 1998; 
Baraffe et al. 2002). Although the BT-Settl models are com- 
puted using more complete line lists than previous genera- 
tions of phoenix models, there is still a need for additional 
data to provide fully comprehensive line lists, especially for 
the molecular opacities (H2O, TiO, CH4, etc.; see Partridge 
& Schwenke 1997 and Schwenke 1998). 

If we are to be quantitative about the missing opac- 
ity we need to measure the missing flux individually in 
each bandpass, rather than as a function of colour. The 
most straightforward way of achieving this is to identify a 
photometric bandpass where the missing sources of opacity 
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are minimal, and create colours with respect to that band. 
This is equivalent to using the luminosity in that band to 
define the temperature. To test the temperature scale re- 
quires masses (from which the models predict temperatures) 
from binaries, data which are simply not available for the 
Pleiades. Instead, as the opacity problem affects MS stars 
as well as pre-MS stars (in the same T e ff range), we will use 
MS binaries to identify a band we can use. 



5 COMPARING THE MODELS AND THE 
DATA - MAIN-SEQUENCE BINARIES 

Although there are many MS binaries with well determined 
masses and distances we could use to test the models, the 
outstanding problem is the small sub-sample of these which 
have individual magnitudes in many colours. Astrometric 
binaries are often separated in the infrared (IR) , and so lack 
individual optical colours, and eclipsing binaries often have 
light curves in only a small number of colours. To overcome 
this problem we will model the system magnitude, which 
as we shall discuss later still gives significant insight as a 
function of stellar temperature. 



5.1 The sample 

To create our sample, data for low-mass binaries from 
Delfosse et al. (2000) were supplemented with data for 
higher-mass binaries from Andersen et al. (1987); Clausen 
et al. (2008); Clausen et al. (2009); Clement et al. (1997b); 
Clement et al. (1997a); Lacy et al. (1997); Lacy et al. (2005); 
Lopez-Morales & Ribas (2005); Popper et al. (1986); Pop- 
per (1994); Popper (1997); Torres et al. (1997), and Torres 
et al. (2009). We imposed an upper limit to the mass un- 
certainty on both the primary and secondary components 
of 10 per cent, and in addition limited ourselves to non- 
contact binaries in which the effects of tidal forces and mass 
transfer may bias the results. We also limited the sample 
to binaries with associated uncertainties in the parallax of 
< 11 per cent; thus our conclusions are unaffected by the 
Lutz-Kelker bias (Lutz & Kelker 1973). The sample span 
the mass range 0.1 < Mq < 1.4 for the individual compo- 
nents, with the higher mass limit chosen as the BCAH98 
models only extend to such masses. For the low-mass sys- 
tems in Delfosse et al. (2000), the apparent magnitudes were 
from the homogeneous dataset of Leggett (1992). A number 
of higher mass systems only had a Stromgren b — y colour 
and so these were transformed into V — I c colours according 
to the relation of Bessell (1979) which is accurate to within 
0.01 mag for MS stars. Near-IR measurements were taken 
from the 2MASS point source catalogue (Cutri et al. 2003). 
Parallax measurements were taken from the Hipparcos cata- 
logue (Perryman et al. 1997), the Yale General Catalogue of 
Trigonometric Stellar Parallaxes (van Altena, Lee, & Homeit 
1995) and Segransan et al. (2000). 



5.2 The models 

We adopted 2 Gyr isochrones with which to compare obser- 
vations against theory as we expect these binaries, which 
are located in the field, to be old systems. Note that at such 
ages, the mass scale is insensitive to changes in age. For 
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Figure 10. The difference between the observed and theoretically 
predicted system absolute magnitudes (where A implies calcu- 
lated - observed) for the 28 binary systems as a function of the 
measured primary component stellar mass. Bolometric luminosi- 
ties have been calculated using 2 Gyr interior models of BCAH98 
a = 1.9 (red), SDF00 (blue), DAM97 (green) and DCJ08 (cyan). 
The atmospheric model flux distribution were folded through the 
revised UBVRI responses of Bessell Sz Murphy (2012) and the 
JHK S responses of Cohen et al. (2003) to create bolometric cor- 
rections which were used to derive the predicted absolute magni- 
tudes. Upper panel: AV. Upper middle panel: A/ c - Middle panel: 
A J. Lower middle panel: AH. Lower panel: AK S . 



each binary component, the bolometric luminosity (Lboi), 
defined by the stellar interior models, was converted into an 
absolute magnitude using bolometric corrections calculated 
using Eqn. B2 by folding the atmospheric models through 
the revised UBVRI responses of Bessell & Murphy (2012), 
and the JHK S responses of Cohen, Wheaton, & Megeath 
(2003). For the reference Vega spectrum we used the CAL- 
SPEC alpha_lyr_stis_005 4 spectrum with V = 0.03 and all 
colours equal to zero. For the 2MASS responses we used 
the bandpass specific reference spectra and zero-point off- 
sets of Cohen et al. (2003) (see their Table 2 and Section 5 
respectively). The absolute magnitude was then converted 
into a flux, and the fluxes for a given binary pair were then 
summed and converted back into an absolute system magni- 
tude. Where possible, we have included the effects of inter- 
stellar reddening by adopting literature values. These values 
range from E(B — V) = 0.0 — 0.11. The atmospheric models 



http : //www . st sci . edu/hst/observatory/cdbs/calspec . html 



10 Cameron P. M. Bell et al 



were reddened accordingly and the bolometric corrections 
calculated as discussed in Section 4.1. 

5.3 Discussion 

Fig. 10 shows the residuals (calculated — observed) in lumi- 
nosity relation in different bandpasses. As the binary mass 
ratio q —> 1 both components have the same mass, whereas 
when q — >• the primary component is likely to dominate 
the flux output of the system. Thus we plot the difference 
in absolute magnitude as a function of the measured pri- 
mary mass. The range of q as defined by our sample is 
q — 1712/ mi = 0.33 — 1.0, where for 65 per cent of these 
the ratio q > 0.9. Regardless of the choice of interior model, 
there are significant systematic deviations (especially in the 
optical regime) for masses below ~ 0.4 — 0.5 M©. Fig. 10 
shows that the differences between the observed and calcu- 
lated system absolute magnitude decreases as one moves to 
redder bandpasses, from approximately 2.5 mag in the V- 
band to less than 1 mag in the if s -band. In the near-IR, and 
in particular the i^ s -band, there is far less spread for a given 
stellar mass (see also Martin et al. 2000 and Delfosse et al. 
2000), with the interior models predicting similar luminosi- 
ties. This further highlights the fact that there is a problem 
with the optical colours predicted by the atmospheric mod- 
els for masses below ~ 0.4 - 0.5 M (T eff < 3700 K). 

It has been noted in the literature that short-period, 
close-in binary system components have radii which are in- 
flated of order 10 per cent with respect to evolutionary mod- 
els of low-mass MS stars (e.g. Kraus et al. 2011). To test 
whether we observe this effect in the photometry of these 
systems, we plotted Fig. 10 as a function of physical separa- 
tion. We found no evidence for the models to either under- or 
overpredict the absolute system magnitude for short-period 
binaries with smaller physical separations. 

What is clear, however, is that of all the bands available, 
the ifg-band magnitude is closest to that predicted by the 
models for a given temperature, both in terms of mean mag- 
nitude and spread. This is unsurprising, since the if s -band 
is less affected than either the J- or if-band by the missing 
sources of opacity that arise from incomplete H2O line lists 
(see the discussion on the completeness of the BT2 H2O line 
list as a function of wavelength by Barber et al. 2006). In 
what follows, therefore, we will use the i^ s -band magnitude 
to determine the temperatures of Pleiades members. 



6 QUANTIFYING THE DISCREPANCY 

Our approach for quantifying the effects of the missing opac- 
ity is to determine the temperature of a given Pleiad by 
comparing its i^-band magnitude with that predicted by a 
given model. We can then compare the predicted luminosity 
in any other band with that of the member in question to 
calculate the missing flux. In practice it is better to work 
with a sequence than individual members and so we de- 
fined the single-star Pleiades sequence in the INT-WFC and 
2MASS colours by fitting a spline (by eye) to the observed 
sequence in various CMDs. We used only stars with uncer- 
tainties in both colour and magnitude of less than 0.1 mag. 
This spline lies slightly above the lower envelope of the se- 
quence to account for photometric uncertainties and takes 



into consideration that the equal-mass binary sequence lies 
~ 0.75 mag above the single-star sequence (with very few 
systems of higher multiplicity) . Table 6 provides the single- 
star sequence in the (griZ) WFC JHK s bandpasses. 

Starting, for instance, with the K s , (g WFC —K s ) Pleiades 
CMD, we compared the theoretical 130 Myr isochrone to 
the observed spline and thus calculated the difference be- 
tween the colours at a given T e ff (provided by the model 
isochrone). Assuming that the problem, in this case, lies in 
the g WFC -band and not K s , this then defines the required 
correction (ABC; defined as BCobs — BC ca ic) for the g WFC - 
band at a specific T e ff. This process was then repeated over 
the T e ff range of the theoretical isochrone or the colour range 
specified by the observed spline, whichever was more re- 
strictive, thus providing ABC as a function of T e ff for the 
g WFC -band. Repeating this for all four pre-MS evolutionary 
models (including both the BCAH98 computations), we cal- 
culated the model dependent ABCs as a function of T e ff for 
the (griZ) WFC JH bandpasses. 

Our final results are shown in Fig. 11. In the Pleiades, 
stars cooler than T e ff ~ 4000 K are on the pre-MS, and all 
models fail to match the data in the optical in this regime. 
The discrepancy is large, 0.75 mag in g WFC , but decreases 
with increasing wavelength. This is a clear improvement over 
the PHOENix/GAIA models, for which we find a discrepancy 
of around 1 mag. 

We can test what effect the overestimation in flux of 
the atmospheric models has on the derived ages for T e ff < 
4000 K in the following way. If this discrepancy is indepen- 
dent of age (and therefore logg), the roughly logarithmic 
age spacing between the isochrones in CMD space means 
that the resulting error in age is best represented as a frac- 
tional rather than an absolute difference. We begin by tak- 
ing the BCAH98 a = 1.9 isochrones and adding our ABCs 
derived above to the theoretical bolometric corrections and 
colour- T e ff relations at a given age. We then find the closest 
matching uncorrected BCAH98 a = 1.9 isochrone (in CMD 
space). Comparing these, we find that for ages between 3 
and 10 Myr, the difference in age is approximately a factor 
3 in g, g — i and 2 in r, r — z. 

The message here is very clear, one should use the red- 
dest waveband possible to determine an age. Unfortunately 
by the time one reaches the JHK S bandpasses, the pre- 
MS isochrones are vertical with a colour J — K s c± 0.85 
for 2500 < Teff < 4000 K. Although individual stars descend 
with age, the resulting sequences are almost degenerate with 
age, and for young clusters (< 10 Myr), observations are fur- 
ther complicated by the presence of discs. Thus the JHK S 
IR is unlikely to yield reliable young cluster ages. Masses 
do not suffer from such a degeneracy and therefore the H- 
band may be a reliable mass indicator when K s is affected 
by circumstellar material. 

Perhaps more surprising is the fact that in the MS 
regime, hotter than T e ff ~ 4000 K, only two models fit the 
data in the optical. The reason the BCAH98 a = 1.0 and 
SDF00 models tend to overpredict the luminosity at a given 
T e ff is that neither is tuned to the Sun, as explained in Sec- 
tion 2.1. Once the BCAH98 models have a tuned mixing 
length paramater (a = 1.9) based on matching observed so- 
lar values they are a good match to the MS data. A possible 
reason why the DAM97 models fail to fit the MS may be 
because of their treatment of convection. 
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Figure 11. Model dependent corrections (ABC) calculated as a function of T e ff for the optical (griZ) WFC and near-IR JH bandpasses. 
The corrections were calculated as the difference between the theoretically predicted colour and the observed colour of the Pleiades 
sequence at a given T e ff. Corrections have been calculated for the following models: BCAH98 a = 1.9 (red, continuous), BCAH98 
a = 1.0 (red, dashed), SDF00 (blue), DAM97 (green) and DCJ08 (cyan). For T eff < 4000 K all the models fail to match the data in the 
optical regime, with the magnitude of the mismatch decreasing with increasing wavelength. 



Table 6. The Pleiades single-star sequence in the INT-WFC and 
2MASS bandpasses. The full table is available online, a sample is 
shown here as a representation of its content. 



^WFC 


r WFC 


Vfc 


^WFC 


J 


H 


K s 


10.807 


10.377 


10.270 


10.251 


9.370 


9.082 


8.987 


10.972 


10.505 


10.389 


10.364 


9.484 


9.181 


9.078 


11.133 


10.645 


10.516 


10.481 


9.594 


9.275 


9.166 


11.292 


10.783 


10.640 


10.596 


9.699 


9.364 


9.251 


11.445 


10.914 


10.758 


10.704 


9.796 


9.445 


9.330 


11.593 


11.040 


10.871 


10.809 


9.887 


9.520 


9.405 


11.735 


11.161 


10.979 


10.908 


9.970 


9.588 


9.473 


11.870 


11.273 


11.079 


11.000 


10.043 


9.646 


9.534 


11.998 


11.378 


11.172 


11.086 


10.110 


9.699 


9.589 


12.120 


11.478 


11.261 


11.168 


10.169 


9.745 


9.637 



(i) We have demonstrated that traditional photometric 
calibration, using observations of standard stars to trans- 
form the data into a standard photometric system, should 
not be used for studies of pre-MS stars. Differences in logg 
and T e ff between a MS and pre-MS star of the same colour 
can, for instance, result in errors of order 0.1 — 0.15 mag 
in both the g-band magnitude and g — i colour for red 
stars. Hence it is crucial that precise photometric studies 
(especially of pre-MS objects) be carried out in the natural 
photometric system of the observations. 



(ii) Therefore we have calculated system responses for the 
natural INT-WFC photometric system. Our observations 
show these are a good model of the photometric system. 



7 CONCLUSIONS 

We have carried out a precise test of a set of pre-MS models 
by examining how well they match the Pleiades sequence in 
optical and near-IR CMDs. The stages we have gone through 
to achieve this are as follows. 



(iii) We have demonstrated that for all optical colours no 
pre-MS model follows the observed Pleiades sequence for 
temperatures cooler than 4000 K. The models overestimate 
the flux by a factor two at 0.5 /xm, with the difference 
decreasing with increasing wavelength. We believe there is 
little observable effect at 2.2 /xm, as the if s -band magnitude 
of MS binaries matches the models at this wavelength. 
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(iv) These differences in bolometric correction corre- 
spond to underestimating the ages of pre-MS stars younger 
than lOMyr by factors of two to three, depending on the 
choice of colours used. Thus the errors in pre-MS isochrone 
models can explain the discrepancy between MS and 
pre-MS ages found by Naylor (2009), and implies the ages 
determined from the MS are more likely to be correct. 

(v) We have provided both our system responses and our 
Pleiades sequence in such a way that it can be used as 
a benchmark for any new evolutionary models. These are 
given in Tables Al and 6 respectively and will also be made 
available from the Cluster Collaboration homepage 5 and the 
CDS archive. 
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APPENDIX A: THE CALCULATED INT-WFC 
SYSTEM RESPONSES 

Although the SDSS survey system responses are well con- 
strained and available on the SDSS webpages 6 (Doi et al. 
2010), the responses on the Isaac Newton Group (ING) web- 
pages for the WFC only combine the filter throughput and 
CCD quantum efficiency. To model the system responses 
of the INT-WFC we included the cumulative effects of the 
transmission of the Earth's atmosphere, the reflectivity of 
the telescope mirror, the transmission of the prime focus 
corrector optics, the quantum efficiency of the detector, and 
the filter transmission. To calculate the atmospheric trans- 
mission we used the model for the La Palma atmosphere 
derived by King (1985), for an airmass typical for our ob- 
servations of 1.4. This model varies smoothly as a function 
of wavelength and does not include the molecular absorp- 
tion features. These atmospheric absorption bands (primar- 
ily due to H2O, CO2 and O3) were estimated using the 
spectrum of an F8 star observed using the Faint Object 
Spectrograph on the William Herschel Telescope. At low 
resolution the continuum of an F8 star is relatively smooth 
and thus the atmospheric bands can be identified by fitting 
a low order polynomial to the continuum and subtracting 
this from the spectrum (Shahbaz et al. 1996). The contin- 
uum is modelled between 4700 — 9800 A. At prime focus the 
INT optical path involves a single reflection from the alu- 
minium coated primary mirror and this was modelled us- 
ing the aluminium reflectivity spectrum of Allen (1963). A 
three-element prime focus corrector is used which is coated 
to minimise reflection and thus improve efficiency and re- 
duce ghosts. The first two elements are no n- interchangeable 
and have a broadband single layer MgF2 coating. The third 
element, although interchangeable, typically adds an addi- 
tional single layer of MgF 2 coating. These single layer coat- 
ings produce total reflectivities smaller than 2 per cent in 
the wavelength region 3500 — 7400 A. Beyond this range, 
the reflectivity from the prime focus corrector increases to 
c± 3 per cent at 8000 A. The filter and detector responses 
were from the ING webpages 7 . The detector is comprised 
of four EEV42-80 CCDs. The detector response on the ING 
webpages only extends to 3400 A in the blue, however the 
blue edge of the RGO U WFC filter is at approximately 3000 A 
and so the detector response was extended blueward using 
data for an almost identical EEV44-82 CCD to a wavelength 
of 3200 A and further extrapolated to 3000 A (Cavadore & 
Dorn 2000). The RGO Z WFC filter response data becomes 
noisy at wavelengths greater than 9000 A and so this re- 
sponse was extended redward using calculated 4-mm Schott 

6 http : / / www . sdss . org/ dr7/ instruments/ imager/ index . html# 
filters 

7 http: //www. ing. iac . es/Astronomy/instruments/wf c/ 
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Figure Al. The normalised system response functions for the 
Sloan Digital Sky Survey (dashed) and the Wide-Field Camera 
on the Isaac Newton Telescope (INT-WFC; bold). The INT-WFC 
system responses are those calculated in Appendix A and include 
the effects of the telescope optics and atmospheric absorption. 
Plotted in the bottom panel are normalised model spectra of Vega 
(A0V), the Sun (G2V) and Proxima Centauri (M6V). 



RG850 data to ~ 10 300 A where the detector response ap- 
proaches zero. Any negative values for the filter responses 
from the ING webpages were removed. Note also that the 
RGO U WFC filter response used here does not include the re- 
cently discovered red leak near 7050 A (Verbeek et al. 2012). 
The calculated INT-WFC system responses are shown in the 
top panel of Fig. Al, with the individual bandpass through- 
puts given in Table Al. 



APPENDIX B: CALCULATING THE 
THEORETICAL TRANSFORMATIONS 

To investigate what effect variations in the system responses 
between the standard SDSS and natural INT-WFC photo- 
metric systems would have on the calibration of our pho- 
tometric data, we calculated the magnitude and colour dif- 
ferences between the two photometric systems. We adopted 
the standard SDSS responses, the calculated INT-WFC re- 
sponses and folded the atmospheric models through both. 
For this comparison we were interested only in the relative 
differences in the derived colours and magnitudes, and so 
our conclusions are not be affected by the use of atmospheric 
models. 

Following the formalism of Girardi et al. (2002), a syn- 
thetic apparent magnitude m Rx in a given filter with re- 
sponse function R\ is calculated as, 



= -2.5 log 



f x XfxRxdX 
f x \f° x Rxd\ 



(Bl) 



where fx is the stellar flux as observed at Earth and fx 
denotes a reference spectrum with a known apparent mag- 
nitude m° R . The response function Rx is the product of the 
component responses described in Appendix A. Note that 
the integrands in Eqn. Bl imply photon counting across the 
filter which is appropriate for detectors such as CCDs. 



Stellar atmospheric models only contain information 
about the flux at the stellar surface, Fx. The transforma- 
tion from fx to Fx requires knowledge of the stellar radius, 
its distance and an extinction curve Ax to account for in- 
terstellar extinction. To incorporate these and remove the 
radial dependence, it is more convenient to deal in terms of 
bolometric corrections, which can be formulated as, 



Mboi,© - 2.5 log 



2.5 log 



/ 47r(10pc) 2 Fboi 
f x \F x 10-° AA *Rx dX 



(B2) 



where Fboi = crT eS is the total flux emergent at the stellar 
surface and all other symbols retain their original definitions. 
For the solar values we used Mboi,© = 4.74 and L© = 3.855 x 
lO^ergs -1 (Bessell, Castelli, & Plez 1998). 

The SDSS photometric system is based on a monochro- 
matic AB magnitude system (see Oke & Gunn 1983; 
Fukugita et al. 1996) which by definition means that a ref- 
erence spectrum of constant flux density per unit frequency 
f° will have AB magnitudes ra° = at all frequencies v. By 
converting this reference spectrum in Eqn. B2 to energy per 
unit wavelength according to fx — c/\ 2 f„, this definition 
can then be extended to any photometric system with the 
result that m° x = are the reference magnitudes. Thus the 
magnitude and colour differences for a given flux distribu- 
tion between the INT-WFC and SDSS photometric systems 
were then simply calculated as the difference in the derived 
bolometric corrections. 

We required bolometric corrections for the MS and pre- 
MS regimes. To model the changing T e ff and logg values 
along the MS (from logg c± 4 for early B-type stars to c± 
5.5 for late M-type) we used the calculated ZAMS of SDF00 
for stellar masses in the range 0.1 — 7 Mq. For the pre-MS 
values we adopted the SDF00 interior models for ages of 1 
and lOMyr, calculating T e ff and logg over the same mass 
range. The reason for choosing the SDF00 interior models 
over other pre-MS stellar interior models was simply because 
they span a greater range of stellar masses and so more 
thoroughly map the MS. Due to the degree of model depen- 
dency in pre-MS evolutionary models, we repeated the cal- 
culation of the theoretical transformations using the ZAMS 
Teff, \ogg values from the BCAH98, DAM97 and DCJ08 
models and found that all four sets of evolutionary models 
agree to within 0.02 mag over the range of stellar masses. 

Fig. 3 shows the calculated magnitude and colour differ- 
ences for reddened and unreddened objects (for illustrative 
purposes we adopted a nominal reddening of E(B — V) = 1). 
The simple formalism adopted in the derivation of the bolo- 
metric corrections allowed us to include the effects of in- 
terstellar reddening by applying a given extinction curve 
Ax and calculating the bolometric correction as shown in 
Eqn. B2. We adopted the Cardelli et al. (1989) extinction 
curve of Rv = Ay / 'E(B — V) and thus derived the extinc- 
tion in a given bandpass at a specific E(B — V). We adopted 
Rv = 3.2 as this value best reproduced the colour dependent 
reddening vector of Bessell et al. (1998) upon folding the 
ATLAS 9 "no-overshoot" atmospheric models (Kurucz 1992) 
through the UBVRI responses of Bessell (1990). Magnitude 
and colour differences were calculated for each T e ff, logg 
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Table Al. The normalised calculated INT-WFC system responses. 
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combination in both the MS and pre-MS regimes via linear 
interpolation. 



